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In a scalar field theory, when the tree level potential admits broken symmetry 



' ground states, the quantum corrections to the static effective potential are complex. 

+-> I (The imaginary part is a consequence of an instability towards phase separation 

and the static effective potential is not a relevant quantity for understanding the 
1^ ' dynamics). Instead, we study here the equations of motion obtained from the one 

loop effective action for slow rollover out of equilibrium. We considering the case in 
which a scalar field theory undergoes a rapid phase transition from Tj > Tc to Tj < T^. 
We find that, for slow rollover initial conditions (the field near the maximum of the 
tree level potential), the process of phase separation controlled by unstable long- 
wavelength fluctuations introduces dramatic corrections to the dynamical evolution 
of the field. We find that these effects slow the rollover even further thus delaying 
the phase transition, and increasing the time that the field spends near the "false 
vacuum". Moreover, when the initial value of the field is very close to zero, the 
dynamics becomes non-perturbative. 
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I. INTRODUCTION AND MOTIVATION 



Inflationary cosmological models provide very appealing scenarios to describe the early 
evolution of our universe Since the original model proposed by Guth several alter- 
native scenarios have been proposed to overcome some of the difficulties with the original 
proposal. 

Among them, the new inflationary model is perhaps one of the most attractive. The 
essential ingredient in the new inflationary model is a scalar field that undergoes a second 
order phase transition from a high temperature symmetric phase to a low temperature broken 
symmetry phase. The expectation value (or thermal average) of the scalar field (j) serves as 
the order parameter. Initially at high temperatures, the scalar field is assumed to be in 
thermal equilibrium and 0^0. The usual field- theoretic tool to study the phase transition 
is the effective potential 0-0. At high temperatures, the global minimum of the effective 
potential is at = 0, whereas at low temperatures there are two degenerate minima. 

The behavior of the phase transition in the new inflationary model is the following: as the 
universe cools down the expectation value of the scalar field remains close to zero until the 
temperature becomes smaller than the critical temperature. Below the critical temperature, 
when the effective potential develops degenerate minima away from the origin, the scalar 
field begins to "roll down the potential hill". In the new inflationary scenario, the effective 
potential below the critical temperature is extremely flat near the maximum, and the scalar 
field remains near the origin, i.e. the false vacuum for a very long time and eventually rolls 
down the hill very slowly. This scenario thus receives the name of "slow rollover". During 
the slow rollover stage, the energy density of the universe is dominated by the constant 
vacuum energy density = 0), and the universe evolves rapidly into a de Sitter space 



(see for example the reviews by Kolb and Turner |10], Linde [11] and Brandenberger WW )- 
Perhaps the most remarkable consequence of the new inflationary scenario and the slow 
rollover transition is that they provide a calculational framework for the prediction of density 



fluctuations [jT3[. The coupling constant in the typical zero temperature potentials must be 
fine tuned to a very small value to reproduce the observed limits on density fluctuations 

SHI- 



This picture of the slow rollover scenario, is based on the static effective potential. The 
use of the static effective potential to describe a time dependent situation has been criticized 



by Mazenko, Unruh and Wald [|14|. These authors argued that the dynamics of the cooling 
down process is very similar to the process of phase separation in statistical mechanics. They 
argued that the system will form domains and that the scalar field will relax to the values 
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at the minima of the potential very quickly. 

Guth and Pi |T5[, performed a thorough analysis of the effects of quantum fluctuations 
on the time evolution. These authors analyzed the situation below the critical temperature 
by treating the potential near the origin as an inverted harmonic oscillator. They recognized 
that the instabilities associated with these upside-down oscillators lead to an exponential 
growth of the quantum fluctuations at long times and to a classical description of the proba- 
bility distribution function. Guth and Pi also recognized that the static effective potential is 
not appropriate to describe the dynamics, that must be treated as a time dependent process. 

Subsequently, Weinberg and Wu [0, have studied the effective potential, particularly in 
the situation when the tree level potential allows for broken symmetry ground states. These 
authors carefully analyzed the contributions to the effective potential and showed that the 
imaginary part of the one-loop effective potential is the result of an analytic continuation 
of the unstable modes (inverted oscillators) studied by Guth and Pi. This imaginary part 
is in fact concealing the growth of this unstable modes and a time dependent situation that 
cannot be described by equilibrium statistical mechanics. The imaginary part of the effective 
potential was related to the lifetime of a particular initial quantum state. 

There have been several attempts to study the time evolution of the scalar field either 
in fiat spacetime or de Sitter space |1T7|-^, but to our knowledge, the influence of the 
instabilities that are responsible for domain growth and phase separation, on the dynamics 
of the scalar field has not yet been elucidated. 

In this article we study the quantum dynamics of the scalar field by analyzing the situ- 
ation in which the system originally in equilibrium at a temperature higher that the critical 
temperature evolves out of equilibrium at a temperature below the critical temperature. We 
provide a detailed description of the dynamics out of equilibrium, concentrating on the in- 
stabilities that drive the phase transition and phase separation, that is the unstable growth 
of long-wavelength fluctuations. Our approach consists in obtaining the non-equilibrium 
evolution equation obtained from the one-loop effective action. These equations turn out to 
be non-local and non-integrable [Q. We provide a qualitative discussion and a consistent 
numerical analysis of these evolution equations for a wide variety of initial conditions. 

The results of our analysis show that the instabilities that trigger the growth of long- 
wavelength fluctuations dramatically enhance the quantum corrections and result in a further 
slow-down of the scalar field. 

For "slow-rollover" initial conditions, we show that even for very weak coupling (consis- 
tent with the bounds on density fluctuations), the quantum corrections become very impor- 
tant, and slow the dynamics dramatically, and in particular for some initial conditions, in 
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which the scalar field is very close to the "false vacuum", the quantum corrections must be 
treated beyond perturbation theory. 

This paper is organized as follows: in the next section we review the effective potential 
and its shortcomings to describe the dynamics. 

In section (111) we describe the non-equilibrium formalism and our approach to obtain 
the equations of motion out of equilibrium for the spatial average of the scalar field. In this 
section, we provide an analysis of the renormalization aspects of the effective equations of 
motion. 

In section (IV) we study the dynamics and provide an analytic as well as numerical 
analysis of the evolution equations for a wide range of initial conditions emphasizing the 
consequences of the unstable growth of long wavelength fluctuations. 

We conclude the paper with an analysis of the potential implications of our results on 
cosmological phase transitions and inflationary models. 



II. THE COMPLEX EFFECTIVE POTENTIAL 



It is already well known that when the tree level potential allows for broken symmetry 
ground states, the one loop effective potential becomes complex in a region of fleld conflgura- 
tions. This unacceptable complex part of the effective potential has been deemed an artifact 
of the loop expansion, and for the most part ignored in most treatments of the effective 
potential. 



However as was clearly shown by Weinberg and Wu |T6|, the imaginary part of the 
effective potential has a very physical meaning, and is a consequence of the instabilities that 
drive phase separation. 

A rather clear understanding of the imaginary part and the physics associated with it is 
obtained by a derivation of the effective potential within the Hamiltonian framework. The 
Hamiltonian for a scalar fleld theory quantized in a volume Q is 

H = l^£x{^n'ix) + ^iV^x)Y + Vi^)^ (2.1) 

Since the effective potential is a function of the zero momentum component of the field, 
we separate the constant part (zero momentum) {(f) of the field and its canonical momentum 

$(x) = (^ + ^(x) ; if=l-(d^x^ix)] [d^xij(x) = (2.2) 

il Jn Jn 

U{x) = ^ + n{x) ; P= f d^xn{x) ; / d^xn{x) = (2.3) 
il Jn Jn 

[U{x), <l>{y)] = -ih5\x - y) ; [P, ^p] = -ih (2.4) 
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We have kept (h) in the above commutators to clarify the quantum corrections. 
Using that n jip, do not have zero momentum component, the Hamiltonian becomes 

H = ^ + nViv) + d?x {^vr^(x) + \m{x)f + \v"{^)i;\x) + ■ ■ j (2.5) 

where the elhpsis stand for higher order terms. 

To this order the Hamihonian for the field il){x) is quadratic, with a mass term that de- 
pends on if. In the representation with creation and annihilation operators for the harmonic 
oscillators, the Hamiltonian finally becomes (in the discrete momentum representation) 

H=^ + ^V{v) + ^(4a, + ]^)hu^{k, if) (2.6) 

k 

u;{k,^) = [e + V"{^)]^ (2.7) 

When all the oscillators are in their ground state, the zero point energy for the oscillators is 
then recognized as the first order quantum correction to the effective potential. Taking the 
large volume limit, the effective potential to this order is then 

Veffi^) = V{^) + 2 / + ^"(^)] ' (2-8) 

As usual renormalization is carried out in the standard manner. From now on, we set h = 1. 
Consider for simplicity the case in which 

Vi^) = -^f^l^' + (2.9) 

then for all values of ip such that V"{ip) < 0, the effective potential acquires an imaginary 
part given by 

ImV.ffiif) = /'"'""^ e \\V"{ip)\ - e] ^ dk (2.10) 



An"^ JO 

{kma.f=\V"{^)\ (2.11) 

This imaginary part is finite (independent of renormalization), and arises because for 
these values of ip the modes for which k"^ < {kmaxY are unstable, the frequencies of the 
oscillators in ( |2.6| ) are imaginary. For these modes the potential is that of an inverted 
oscillator. 

These inverted oscillators for which oj'^ik^ip) < 0, do not have stationary state solutions 
and the imaginary part of the zero point energy arises from the analytic continuation of har- 
monic oscillator wavefunctions from positive to negative o;^. The sign of the imaginary part 
depends on the analytic continuation. The quantum mechanics of these inverted harmonic 



oscillators has been thoroughly studied by Guth and Pi [|15| and Weinberg and Wu |]T6| 
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There being no stationary states associated with these inverted oscillators, an initial state 
must be specified. If, for example, a gaussian wave packet is prepared initially centered at 
the origin of these inverted oscillators, it will spread and the width of the packet will increase 



exponentially in time |]15|,|T6[. This exponential growth, is manifest in the equal time two 
point function [p!5| , p!^ 



This is the two-point Green's function computed in this Gaussian state and measures the 
fiuctuations of the operator ip. In particular, the growth of these unstable modes is precisely 
the mechanism that Mazenko et. al. [0] suggest. 

One notes, however, that the two-point Green's function calculated in the analytically 
continued state is time independent and given by 

<^,{t)^_,{t)>=^-^^ (2.12) 

This result is purely imaginary, but again this imaginary part is obtained after the ana- 
lytic continuation of these Gaussian states, and is actually concealing a time dependent 
situation. 

Another crucial observation is that it is this two-point function evaluated in these Gaus- 
sian states that gives the one-loop correction to the effective potential. Thus, it becomes 
clear, that the imaginary part, resulting from an analytic continuation of the unstable modes 
is in fact hiding an unstable time evolution. At finite temperature, the imaginary part is 
signaling a non- equilibrium situation. 

This situation is well understood in statistical mechanics, perhaps the earliest example is 
the Van der Walls equation of state and its unphysical isotherms. In the unphysical region, 
the situation must be studied out of equilibrium. In this region there is coexistence of 
different phases that may not be studied within equilibrium statistical mechanics. An ad-hoc 
remedy in this situation, is the Maxwell construction that replaces the unphysical region of 
the isotherms by a straight line. Its physical interpretation is that the non-equilibrium state 
may be found, in the coexistence region, as a mixture of phases with arbitrary concentrations 
of each phase. This situation is also typical of binary mixtures in statistical mechanics 
p3| , |2^ . The Maxwell constructed effective potential or free energy is irrelevant for the 
dynamical non-equilibrium description of the system. 

The exponential growth of the unstable modes (inverted oscillators) is signaling the 
growth of domains and the onset of the process that triggers the phase transition, i.e. 
phase separation. This is very similar to the process of spinodal decomposition in statistical 
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mechanics p3| , p^ . An attempt to describe spinodal decomposition within the context of 
field theory has been described by Calzetta . 

Clearly the approximation of inverted oscillators, is crude as it neglects non-linear effects. 
The growth of these unstable modes will eventually slow down when the non-linearities be- 
come important, this is the process of coarsening. As we will show later, a clear understand- 
ing of the physics of coarsening in the regime where the non-linearities become important, 
may require departing from perturbative treatments. 

Usually in order to understand the dynamics of the scalar field in this type of situations, 
the static effective potential is used in the equations of motion, resulting in a typical evolution 
equation (in fiat spacetime) 

^ + ^%M = o (2.13) 
dt'^ dip 

After the arguments presented above, it becomes clear that this equation is inappropriate 
for the study of the dynamical evolution. The imaginary part of Veff{(p) is signaling an 
unstability and the static effective potential is not a suitable quantity to study the dynamics. 
That is one must consider the effective action for time dependent fields and not merely the 
effective potential which holds for constant fields. 

We turn to this study in the next section. 



III. NON-EQUILIBRIUM TIME EVOLUTION 

As explained in the previous section, one must depart from the usual description in terms 
of the static effective potential, and treat the dynamics with the full time evolution. The 
time evolution of the system will be determined once the evolution Hamiltonian and the 
initial state are prescribed. 

In order to understand the dynamics of the phase transition and the physics of the 
instabilities mentioned above, let us consider the situation in which for time t < the 
system is in equilibrium at an initial temperature Tj > Tc where Tc is the critical temperature. 
At time t = the system is rapidly "quenched" to a final temperature below the critical 
temperature Tf < Tc and evolves thereafter out of equilibrium. 

What we have in mind in this situation, is a cosmological scenario with a period of rapid 
infiation in which the temperature drops very fast compared to typical relaxation times of 
the scalar field. In particular this situation should correspond to the case ^ r with 
h Bubble's constant and r a typical relaxation time. At high temperatures and weakly 
coupled theories we would expect r^^ ^ A^T (A is the coupling). When T ^ Tc oc A~^/^ 
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relaxation times become very large and the dynamics of the long-wavelength modes (the only 
relevant modes for the phase transition) becomes critically slowed down. Precisely because 
of this critical slowing down, we conjecture that an inflationary period at temperatures 
near the critical temperature, may be described in this "quenched" approximation. Another 
situation that may be described by this approximation is that of a scalar field again at 
Tj > Tc suddenly coupled to a "heat bath" at a much lower temperature (below the transition 
temperature) and evolving out of equilibrium. The "heat bath" may be other fields at a 
different temperature. The influence of a "heat bath" in an inflationary universe, has been 
studied in the linearized approximation by Cornwall and Bruinsma . 

One then would expect that a "quenching out of equilibrium scenario" may be an appro- 
priate description near the critical temperature for these situations. Certainly this is only a 
plausibility argument, a deeper understanding of the initial conditions must be pursued to 
obtain a more precise knowledge of the cooling down process. 

We do not envisage in this article to study the case of an inflationary cosmology or 
the detailed dynamics of the mechanism that produces the "quenching" below the critical 
temperature, and departure from equilibrium. We expect to report on these investigations 
in forthcoming articles. 

Here we just assume that such a mechanism takes place and simplify the situation by 
introducing a Hamiltonian with a time dependent mass term to describe this situation 

H{t) = d\ \^-Il\x) + \{V^{x)f + \m\t)^\x) + A$4(^)| (31) 

m2(t) = m'ei-t) + ( V)0W (3.2) 

where both m? and /i^ are positive. We assume that for all times t < there is thermal 
equilibrium at temperature Tj, and the system is described by the density matrix 

p^ = e-^"^' (3.3) 
Hi = H{t < 0) (3.4) 

In the Schroedinger picture, the density matrix evolves in time as 

p{t) = U{t)p,U-\t) (3.5) 

with U (t) the time evolution operator. 

An alternative and equally valid interpretation (and the one that we like best) is that the 
initial condition being considered here is that of a system in equilibrium in the symmetric 
phase, and evolved in time with a Hamiltonian that allows for broken symmetry ground 
states, i.e. the Hamiltonian (ETTl, R.2|) for t > 0. 
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The expectation value of any operator is thus 

<0>{t)= Tre-^^"^U~\t)OU{t) /Tre-^^^^ (3.6) 

This expression may be written in a more illuminating form by choosing an arbitrary time 
T < for which U{T) = exp[—iTHi\ then we may write exp[— = exp[—iHi{T — if3i — 
T)] = U(T — iPi,T). Inserting in the trace U~^(T)U(T) = 1, commuting U^^{T) with pi 
and using the composition property of the evolution operator, we may write ( p.6[ ) as 

<0>{t)= TrUiT - i(3i, t)OU{t, T)/TrU{T - t(3i, T) (3.7) 

The numerator of the above expression has a simple meaning: start at time T < 0, evolve 
to time t, insert the operator O and evolve backwards in time from t to T < 0, and along 
the negative imaginary axis from T to T — ijSi. This operation is depicted in figure 1 (a). 
The denominator, just evolves along the negative imaginary axis from T to T — i(3i. The 
contour in the numerator may be extended to an arbitrary large positive time T' by inserting 
U{t,T')U{T',t) = 1 to the left of O in (g^), thus becoming 

<0>{t)= TrUiT - if3i, T)U{T, T')U{T\ t)OU{t, T)/TrU{T - T) (3.8) 

The numerator now represents the process of evolving from T < to t, inserting the operator 
O, evolving further to T', and backwards from T' to T and down the negative imaginary 
axis to T — ijSi. This process is depicted in the contour of figure 1(b). Eventually we take 
T — oo ; T' —>■ oo. It is straightforward to generalize to real time correlation functions of 
Heisenberg picture operators. 

This formalism allows us also to study the general case in which both the mass and the 
coupling depend on time, and furthermore, by taking the zero temperature limit, we can 
study the situation in which a particular state is prepared at time t = and evolved in time. 

For example, by switching off the coupling for t < one is preparing a gaussian density 
matrix at t = or in the zero temperature limit, a gaussian wavefunctional. This density 
matrix or gaussian functional will then be evolved in time and in this time evolved state (or 
density matrix) we compute expectation values of operators, or correlation functions. We 
then see that the above formalism permits us to study these situations in great generality. 

As mentioned before, another point of view that one may take on the "quenching" 
below the critical temperature, is that a definite state or density matrix describing the 
symmetric phase is prepared as an initial condition for t < and evolved in time with the 
Hamiltonian that allows for broken symmetry states. One then studies the dynamics of 



9 



the phase transition, and how the system evolves in time from the initially symmetric state 
towards the asymmetric states. 

As usual, the insertion of an operator may be achieved by inserting sources in the time 
evolution operators, defining the generating functionals and eventually taking functional 
derivatives with respect to these sources. Notice that we have three evolution operators, 
from T to T', from T', back to T (inverse operator) and from T to T — since each of 
these operators has interactions and we want to use perturbation theory and generate the 
diagrammatics from the generating functionals, we use three different sources. A source 
for the evolution T T', J~ for the branch T' —* T and finally for T ^ T — The 
denominator may be obtained from the numerator by setting J+ = = 0. Finally the 
generating functional Z[J+ , J-, J^] = TrU{T - T; J^)U{T, T'; J-)U{T', T; J+), may be 
written in term of path integrals as (here we neglect the spatial arguments to avoid cluttering 
of notation) 

e^ir''^^^''^'^ (3.9) 

with the boundary conditions $+(T) = $^(T-i/3,) = $ ; $+(T') = <I>-(T') = $2 ; $"(T) = 
$^(T) = $1. As usual the path integrals over the quadratic forms may be done and one 
obtains the final result for the partition function 

els/e'^*! Jjt2Mtl)Mt2)G4h,t2)} (3.10) 



Where Jc are the currents defined on the contour of figure (l,b) , p5[ and Gc is the 



Green's function on the contour (see below), and again the spatial arguments have been 
suppressed. 



In the two contour integrals (on ti, ^2) in (|3.10| ) there are altogether nine terms, cor- 
responding to the combination of currents in each of the three branches. However, in the 
limit T — i> — 00, the contributions arising from the terms in which one current is on the (+) 
or (— ) branch and another on the imaginary time segment (from T to T — i/?j), go to zero 
when computing correlation functions in which the external legs are at finite real time. For 
this real time correlation functions there is no contribution from the terms, that cancel 
between numerator and denominator, and the information on finite temperature is encoded 
in the boundary conditions on the Green's functions (see below). Then for the calculation 



of finite real time correlation functions the generating functional simplifies to [27,28 
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Z[J+, j- 



It '^^l /t' '^*2Ja(tl)Jb(i2)G'„5(ti,t2)} 



X 



(3.11) 



with a,b = +, — . 

This formulation in terms of time evolution along a contour in complex time has been 
used many times in non-equilibrium statistical mechanics. To our knowledge the first to use 
this formulation were Schwinger [^] and Keldysh |30] (see also Mills ||31[]). There are many 
articles in the literature using these techniques to study time dependent problems, some 
of the more clear articles are by Jordan |Q, Niemi and SemenofF Landsman and van 
Weert ||3^,Semenoff and Weiss |Q, Kobes and Kowalski Calzetta and Hu [Q,Paz ||3^ 
and references therein. 

At first sight one seems to have complicated the situation enormously by doubling the 
number of fields. However, this doubling is a natural consequence of dealing with a time 
evolution of a density matrix and in general with probabilities, instead of amplitudes. Rather 
than computing in-out amplitudes, we are here computing expectation values or correlation 
functions in the time evolved in-state or density matrix p9| , p^ . 

The Green's functions that enter in the integrals along the contours in p.lU] , |3.11D are 
given by (see above references) 

G++(tl,t2) = G>{ti,t2)&{t,-t2)+G<{ti,t2)e{t2-t,) (3.12) 
G>{t^, t2)Q{t2 - ti) + G<(ti, t2)e(ti - t2) (3.13) 
-G<{t,M) (3.14) 



G-{t,M) 
G+~{t,M] 
G-^{tiM) 
G<{TM) 



-G^{ti,t2 



-G'^{t2, ti 



(3.15) 

G>{T-iPiM) (3.16) 

As usual G'^^G^ are homogeneous solutions of the quadratic form with appropriate 
boundary conditions. We will construct them explicitly later. The condition ( 3.16 ) is rec- 
ognized as the periodicity condition in imaginary time (KMS condition) It is straight- 
forward to show, using the above Green's functions, that Z[J, J] = 1 as it must. 

Although most of the details presented above on the non- equilibrium formalism are 
available in the literature, we included them here for self-consistency and with the intention 
to clarify some issues that are usually glossed over in most treatments. 

We are now in condition to obtain the evolution equations for the average of the scalar 
field in the case when the potential is suddenly changed, to account for a sudden change 
in temperature from above to below the critical temperature as described by the model 
Hamiltonian (|3.1|) with (|3.2|). For this purpose we use the tadpole method [0], and write 
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(3.17) 



Where, again, the ± refer to the branches for forward and backward time propagation. The 
reason for shifting both (±) fields by the same classical configuration, is that enters in 
the time evolution operator as a background c-number variable, and time evolution forward 
and backwards are now considered in this background. 

The evolution equations are obtained with the tadpole method by expanding the La- 
grangian around (f){t) and considering the linear, cubic, quartic, and higher order terms in 

as perturbations and requiring that 

< ^^(f,t) >= 0. 

It is a straightforward exercise to see that this is equivalent to extremizing the one-loop 
effective action in which the determinant (in the logdet) incorporates the boundary condition 
of equilibrium at time t < at the initial temperature. 

To one loop we find the equation of motion 



{2n) 



-t)Gk{t,t) = 



(3.18) 



with Gkit, t) = Gf(t, t) = G^(t, t) is the spatial Fourier transform of the equal-time Green's 
function. 

At this point, we would like to remind the reader that 

{-tG,{t,t)) =< ^im^_^^{t)> 

is a positive definite quantity (because the field is real) and as we argued before (and 
will be seen explicitly shortly) this Green's function grows in time because of the instabilities 



associated with the phase transition and domain growth |[T5|JT6|1 . 

These Green's functions are constructed out of the homogeneous solutions to the operator 
of quadratic fluctuations 



dt 



^+k^ + M\t) 



Ul 







M\t) = + ^0D0(-^) +( V + ^m)m 

The boundary conditions on the homogeneous solutions are 



Wit < 0) = e 



X 



(3.19) 
(3.20) 

(3.21) 
(3.22) 
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where 0i is the value of the classical field at time t < and is the initial boundary condition 
on the equation of motion. Truly speaking, starting in a fully symmetric phase will force 
(pi = 0, and the time evolution will maintain this value, therefore we admit a small explicit 
symmetry breaking field in the initial density matrix to allow for a small (pi. The introduction 
of this initial condition seems artificial since we are studying the situation of cooling down 
from the symmetric phase. 

However, we recognize that the phase transition from the symmetric phase occurs via 
formation of domains (in the case of a discrete symmetry) inside which the order parameter 
acquires non-zero values. The domains will have the same probability for either value of 
the field and the volume average of the field will remain zero. These domains will grow 
in time, this is the phenomenon of phase separation and spinodal decomposition familiar 
in condensed matter physics. Our evolution equations presumably will apply to the coarse 
grained average of the scalar field inside each of these domains. This average will only 
depend on time. Thus, we interpret ipi as corresponding to the coarse grained average of 
the field in each of these domains. The question of initial conditions on the scalar field is 
also present (but usually overlooked) in the slow-rollover scenarios but as we will see later, 
it plays a fundamental role in the description of the evolution. 

The identification of the initial value (pi with the average of the field in each domain is 
certainly a plausibility argument to justify an initially small asymmetry in the scalar field 
which is necesary for the further evolution of the field, and is consistent with the usual 
assumption within the slow rollover scenario. 

We are currently studying the dynamics of the phase transition from the symmetric 
phase by looking at the composite operator $^(x, t) which measures the fiuctuations, and 
will report on our studies in a forthcoming article . An alternative approach using this 
composite operator has also been proposed by Lawrie [^. 

The boundary conditions on the mode functions U^{t) correspond to "vacuum" boundary 
conditions of positive and negative frequency modes (particles and antiparticles) for t < 0. 

Finite temperature enters through the periodicity conditions ( p.l6| ) and the Green's 
functions are 

Glit^t') = 2u;lik)l-e-^^-<('^) K^^)^^"^^') + e-^'^^^'^U^itMit')] (3-23) 
G<{t,t')=G>{t',t) (3.24) 

Summarizing, the effective equations of motion to one loop that determine the time 
evolution of the scalar field are 
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± 



4>{t) 



(27r)3 2cJ<(A;) 



coth 



(3iUJ<{k) 



(3.25) 
(3.26) 



with (13:201) , (|3:21|) . 

This set of equations is too comphcated to attempt an analytic solution, we will study 
this system numerically shortly. 

However, before studying numerically these equations, one recognizes that there are 
several features of this set of equations that reveal the basic physical aspects of the dynamics 
of the scalar field. 

i) : The effective evolution equations are real. The mode functions V(^{t) are complex 
conjugate of each other as may be seen from the time reversal symmetry of the equations, 
and the boundary conditions ( |3.21D . This situation must be contrasted with the expression 
for the effective potential for the analytically continued modes. 

ii) : Consider the situation in which the initial configuration of the classical field is near 
the origin (pi ^ 0, for t > 0, the modes for which P < {kmaxY {kmaxY = /^^ — | 
unstable. 

In particular, for early times (t > 0), when (pi ^ 0, these unstable modes behave approx- 
imately as 



are 



A, 



1 - I 



W, 



B. = Al 



X 



-k' 



Then the early time behavior of {—iGk{t,t)) is given by 

1 



[-iGkit,t)) 



2uj<{k) 



1 + 



cosh{2Wkt) - 1] 



(3.27) 
(3.28) 

(3.29) 
(3.30) 



coth[piUj<{k)/2] (3.31) 



This early time behavior coincides with the Green's function of Guth and Pi [T^ and Wein 



berg and Wu [|T6l for the inverted harmonic oscillators when our initial state (density matrix) 
is taken into account. 

Our evolution equations, however, permit us to go beyond the early time behavior and 
to incorporate the non-linearities that will eventually shut off the instabilities. 
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These early-stage instabilities and subsequent growth of fluctuations and correlations, are 
the hallmark of the process of phase separation, and precisely the instabilities that trigger 
the phase transition. 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at very early times. 

At early times, the stable modes for which k"^ > {k^axY obtained from ( p. 27]) , (p.28|) 
, ( p.29| ) by the analytic continuation Wk —iuj>{k) = P — /^^ + 



For t < 0, {t)Uk = 1 ^'^d one obtains the usual result for the evolution equation 

dt"^ d(f) 

with Vefficj)) the flnite temperature effective potential but for t < there are no unstable 
modes. 

It becomes clear, however, that for t > there are no static solutions to the evolution 
equations for (pit) ^ 0. 

iii) Coarsening: as the classical expectation value "rolls down" the potential hill, 
increases and {kmax{t)Y = fj? — decreases, and only the very long-wavelength 

modes remain unstable, until for a particular time tg ; {kmax{ts)Y = 0. This occurs when 
4>'^{ts) = 2/i^/A, this is the inflexion point of the tree level potential. In statistical mechanics 



this point is known as the "classical spinodal point" and tg as the "spinodal time" p3|,|2l 



When the classical fleld reaches the spinodal point, all instabilities shut- off. From this point 
on, the dynamics is oscillatory and this period is identifled with the "reheating" stage in 
cosmological scenarios |TT|JT^ . 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at small positive times, as eventually the unstable growth will 
shut-off. 

The value of the spinodal time depends on the initial conditions of (f){t). If the initial 
value (pi is very near the classical spinodal point, t^ will be relatively small and there will 
not be enough time for the unstable modes to grow too much. In this case, the one-loop 
corrections for small coupling constant will remain perturbatively small. On the other hand, 
however, if 0j ^ 0, and the initial velocity is small, it will take a very long time to reach 
the classical spinodal point. In this case the unstable modes may grow dramatically 
making the one-loop corrections non-negligible even for small coupling. These 
initial conditions of small initial fleld and velocity are precisely the "slow rollover" conditions 
that are of interest in cosmological scenarios of "new inflation" . 

Renormalization: 
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As argued above, for t > 0, there are no static solutions to the equation of motion for the 
scalar field. The mode functions U^{t) depend implicitly on the field As in the usual 
situation, one expects ultraviolet divergences in the one-loop correction. It is not clear from 
the equation of motion for the scalar field whether these divergences may be absorbed in a 
redefinition of the mass and coupling constant, or cancelled by local counterterms. Since we 
want to study the time evolution and be able to extract meaningful information, we must 
first understand the renormalization aspects of the effective equation of motion. 

The ultraviolet divergences must be absorbed in m^(t) and A, whose coefficients in the 
equation of motion are and respectively, i.e. time dependent. Since as mentioned 
previously, for t < the situation corresponds to the usual case of the static effective 
potential, renormalization proceeds in the standard manner. However for t > the situation 
is different and to understand it we need the large k behavior of the mode functions. We 
study the short wavelength behavior by a WKB-type analysis: we define e = 1/k and divide 
the equation for the mode functions by k"^ thus obtaining the equation 



+ l + e^M\t) 



(3.32) 



with the boundary conditions (|3.21|) . Let us define the functions V^(t) as the two linearly 
independent solutions for t > to the equation ( p. 3 2| ) with the boundary conditions V^(0"'') = 
1. The mode functions U^{t) solutions to ( p.32| ) with the boundary conditions ( p.21j ) are 
written as 



Ut{t) = c,Vt{t)+dkV,{t) 



(3.33) 
(3.34) 



the coefficients c^, dk are obtained from the matching conditions at t = 0. We propose a 
WKB ansatz for the mode functions V^(t) 



I ■■Jit) 

S±(t) = 5:e"S„^(t) 

n=0 



(3.35) 
(3.36) 



inserting this ansatz in (|3.32|) and comparing powers of e we find the asymptotic behavior 
for large k to be 

V±(t) = e^'(^*+^/o [l - -^(M^it) - M2(0+) 

The leading behavior for large k of the coefficients is found to be 



+ ■ 



(3.37) 
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1 - 



, M2(0+) 
+ 2fc 



M2(0+) 
2k 



+ ■ 



+ 



(3.38) 

(3.39) 
(3.40) 



Inserting these results for the large k behavior in the one-loop contribution, it is straight- 
forward to find that the divergent terms are independent of temperature and we obtain 



(27r)3 2uj^{k) 



coth 



- ^( V + ^0^(t)) In(^) + finite (3.41) 



where A is an upper momentum cutoff, and k a renormalization scale and the finite part is 
time, temperature and k dependent. 

It is clear that these divergences may be cancelled by local counterterms of the usual 
form, where the mass counterterm depends (locally) on time and changes suddenly at t = 0. 
The coupling constant is renormalized in the usual manner. 

It is important to point out that the integral for the one- loop correction may be split 
into the contribution from the unstable modes k"^ < {kmax{t)Y ^^'^ t^iaX: of the stable modes 
k"^ > {kmaxii))"^ ■ It is only the latter that requires renormalization and where the divergences 
reside. The contribution from the unstable modes is finite and does not require renormal- 
izations. The renormalization of the effective action has been done in an alternative manner 



using dimensional renormalization by A van and de Vega |41 



IV. ANALYSIS OF THE EVOLUTION 



As mentioned previously within the context of coarsening, when the initial value of 
the scalar field 0j ~ 0, and the initial temporal derivative is small, the scalar field slowly 
rolls down the potential hill. But during the time while the scalar field remains smaller 
than the "spinodal" value, the unstable modes grow and the one-loop contribution grows 
consequently. For a "slow rollover" condition, the field remains very small (0^(t) <^ 2/i^/A) 
for a long time, and during this time the unstable modes grow exponentially. The stable 
modes, on the other hand give an oscillatory contribution which is bound in time, and for 
weak coupling remains perturbatively small at all times. 

Then for a "slow rollover" situation and for weak coupling, only the unstable modes 
will yield to an important contribution to the one-loop correction. Thus, in the evolution 
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equation for the scalar field, we will keep only the integral over the unstable modes in the 
one loop correction. 

Phenomenologically the coupling constant in these models is bound by the spectrum of 
density fluctuations to be within the range Xr ~ 10^^^ 



10 



-14 



111. The the stable modes 



will always give a perturbative contribution, whereas the unstable modes grow exponentially 
in time thus raising the possibility of giving a non- negligible contribution. 

With the purpose of numerical analysis of the effective equations of motion, it proves 
convenient to introduce the following dimensionless variables 



fj'Rt ; q 



k/ fi 



R 



^R + \^R(t>l 



(4.1) 
(4.2) 



and to account for the change from the initial temperature to the final temperature (Tj > 



Tc ; Tf < Tc) we parametrize PO 



m 



/iij(O) 

^J^R = Aiij(o) 



^ - 1 

c 

1 ^ 

r2 



(4.3) 
(4.4) 



where the subscripts (R) stand for renormalized quantities, and — /i/j(0) is the renormalized 
zero temperature "negative mass squared" and T"^ = 24/i|.(0)/A/j. Furthermore, because 
and Ti > for the unstable modes Ti > {kmax{t)) and we can take 
the high temperature limit coth[/3jC(j<(A;)/2] 2Ti/uj^{k). Finally the effective equations 
of evolution for t > 0, keeping in the one- loop contribution only the unstable modes as 
explained above (g^ < {qmax{'r)Y) become, after using = fi%{q'^ + L"^) 



dr 



^ q^ + L^ ^ 



dr'' 



iqmax(t)) 

{qr. 
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Vox 



R 



27r2 



rp2 
II. 

T2 



(4.5) 

(4.6) 
(4.7) 
(4.8) 



For Ti > Tc and Tf <^ Tc the coupling ( |4.8| ) is bound within the range 10 — 
10~^. The dependence of the coupling with the temperature reflects the fact that at higher 
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temperatures the fluctuations are enhanced. It is then clear, that the contribution from 
the stable modes is always perturbatively small, and only the unstable modes may introduce 
important corrections if they are allowed to grow for a long time. 

From ( [4. 51) we see that the quantum corrections act as a positive dynamical renormal- 
ization of the "negative mass" term that drives the rolling down djTiamics. It is then clear, 
that the quantum corrections tend to slow down the evolution. 

In particular, if the initial value 77(0) is very small, the unstable modes grow for a long 
time before r]{T) reaches the spinodal point rj{Ts) = l/-\/3 at which point the instabilities 
shut off. If this is the case, the quantum corrections introduce substantial modifications to 
the classical equations of motion, thus becoming non- perturbative. If ri{0) is closer to the 
classical spinodal point, the unstable modes do not have time to grow dramatically and the 
quantum corrections are perturbatively small. 

Thus we conclude that the initial conditions on the field determine whether or not the 
quantum corrections are perturbatively small. 

Although the system of equations ( [4.5| , [4.6|) , are coupled, non-linear and integro- 
differential, they may be integrated numerically. Figures (2,3,4) depict the solutions for 
the classical |^ (solid lines) and quantum (dashed lines) evolution (a), the quantum correc- 
tion (b), i.e. the fourth term in ([4.5|) including the coupling g, the classical (solid lines) and 
quantum (dashed lines) velocities (c) and (q'max(T))^ (d). For the numerical integration, 
we have chosen = 1, the results are only weakly dependent on L, and taken g = 10~^, 
we have varied the initial condition r](0) but used ^^^^|^=o = 0. 

We recall, from a previous discussion that //(r) should be identified with the average of 
the field within a domain. We are considering the situation in which this average is very 
small, according to the usual slow-rollover hypothesis, and for which the instabilities are 
stronger. 

In figure (2. a) ri{0) = 2.3 x 10~^ we begin to see that the quantum corrections become 
important at t ~ lO/fiR and slow down the dynamics. By the time that the classical 
evolution reaches the minimum of the classical potential at rj = 1, the quantum evolution 
has just reached the classical spinodal r] = l/VS. We see in figure (2.b) that the quantum 
correction becomes large enough to change the sign of the "mass term" , the field continues 
its evolution towards larger values, however, because the velocity is different from zero, 
attaining a maximum (figure (2.c)) around the time when the quantum correction attains 
its maximum. As 77 gets closer to the classical spinodal point, the unstability shuts off as 
is clearly seen in figure (2.d) and the quantum correction arising from the unstable modes 
become perturbatively small. From the spinodal point onwards, the field evolves towards the 
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minimum and begins to oscillate around it, the quantum correction will be perturbatively 
small, as all the instabilities had shut-off. Higher order corrections, will introduce a damping 
term as quanta may decay in elementary excitations of the true vacuum. 

Figures (3.a-d) show a more marked behavior for t](0) = 2.27 x 10~^ ; ^^^^ = 0, notice 
that the classical evolution of the field has reached beyond the minimum of the potential at 
the time when the quantum evolution has just reached the classical spinodal point. Figure 
(3.b) shows that the quantum correction becomes larger than 1, and dramatically slows 
down the evolution, again because the velocity is different from zero (figure (3.c)) the field 
continues to grow. The velocity reaches a maximum and begins to drop. Once the field 
reaches the spinodal again the instabilities shut-off (figures (3.b,d)) and from this point 
the field will continue to evolve towards the minimum of the potential but the quantum 
corrections will be perturbatively small. 

Figures (4.a-d) show a dramatic behavior for ri{0) = 2.258 x 10"^ ; ^^^^ = 0. The 
unstable modes have enough time to grow so dramatically that the quantum correction 
(figure (4.b)) becomes extremely large ^ 1 (figure(4.b)), overwhelming the "negative mass" 
term near the origin. The dynamical time dependent potential, now becomes a minimum 
at the origin and the quantum evolution begins to oscillate near t] = 0. The contribution of 
the unstable modes has become non-perturbative, and certainly our one- loop approximation 
breaks down. 

As the initial value of the field gets closer to zero, the unstable modes grow for a very 
long time. At this point, we realize, however, that this picture cannot be complete. To see 
this more clearly, consider the case in which the initial state or density matrix corresponds 
exactly to the symmetric case, r/ = is necessarily, by symmetry, a fixed point of the 
equations of motion. Beginning from the symmetric state, the field will always remain at 
the origin and though there will be strong quantum and thermal fluctuations, these are 
symmetric and will sample field configurations with opposite values of the field with equal 
probability. 

In this situation, and according to the picture presented above, one would then expect 
that the unstable modes will grow indefinetely because the scalar field does not roll down 
and will never reach the classical spinodal point thus shutting-off the instabilities. What is 
missing in this picture and the resulting equations of motion is a self-consistent treatment of 
the unstable fluctuations, that must necessarily go beyond one- loop. A more sophisticated 
and clearly non-perturbative scheme must be invoked that will incorporate coarsening, that 
is the shift with time of the unstable modes towards longer wavelength and the eventual 



shutting off of the instabilities. We are currently [38| exploring a Hartree approximation 
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that will incorporate self-consistently these features. Another possible approach would be a 
variational treatment as advocated in references (18-19) or as proposed by Lawrie p9[] . 
Conclusions: 

We have studied the effective equations of motion for the spatially independent average of 
a scalar field evolving out of equilibrium after a second order phase transition. After pointing 
out the severe shortcomings in the usual description in terms of the effective potential, we 
considered the situation in which a scalar field theory is suddenly cooled below the critical 
temperature from an initial state at a temperature higher than the critical temperature. 
We use the tools of statistical mechanics out of equilibrium to study the real time dynamics 
during a "slow rollover stage" . The effective non-equilibrium equations of motion are studied 
both analytically and numerically to one-loop level. We find that the unstable growth of 
long wavelength fluctuations, that trigger the process of phase separation is responsible of a 
very marked behavior in the time evolution of the scalar field. Even for very weak couplings, 
consistent with the bounds from density fluctuations, for the case of "slow rollover" initial 
conditions, the time evolution is dramatically slowed down as a consequence of this unstable 
growth which signal the onset of the phase transition. When the scalar field is very close 
to the origin (at the local maximum of the potential) the unstable modes can grow for 
a long time and the effect of the quantum (and thermal) corrections become very large 
and eventually non-perturbative. We give a qualitative and quantitative description of the 
coarsening process and the eventual shut-off of the instabilities. 

We argue that a comprehensive treatment and understanding of the late stage dynamics 
of the phase transition, and time evolution of the scalar field, involves a non-perturbative 
treatment out of equilibrium. This treatment must incorporate not only the growth but also 
the coarsening features when the initial state is symmetric and the initial value of the scalar 
field is zero, thus remaining zero during the time evolution. 

This is the next step in a consistent analysis of the dynamics, work on this problem is 
currently underway. Clearly, all this must be extended to the case of inflationary cosmology, 
in particular a description of the dynamics out of equilibrium in a de Sitter background is 
the next stage. In this case the physics of coarsening will be complicated by the red-shift of 
wavevectors and more modes enter the unstable regime, clearly a deeper understanding of 
these instabilities is necessary. 

However we believe that the results presented in this article provide a novel insight into 
the long standing problem of the dynamics of a phase transition and the ensuing evolution 
out of equilibrium when quantum and thermal fluctuations are taken into account. It also 
points out that this problem must be studied as a fully time dependent process and that 
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one must abandon the usual treatment in terms of the effective potential. 

The potential consequences for inflationary scenarios are obvious, the instabilities and 
growth of long-wavelength fluctuations enhance the quantum and thermal corrections that 
become large. As a consequence, the "slow-rollover" stage is delayed, the scalar fleld remains 
near the "false vacuum" for a longer period of time giving rise to a longer inflationary stage 
and a delayed completion of the phase transition. 
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Figure Captions: 

Figure 1(a): Contour of evolution in complex time plane, the cross denotes insertion of 
an operator. 

Figure 1(b): Final contour of evolution, eventually T' ^ oo , T — — oo. 

Figure 2(a): rj vs r (notation in the text) for g = 10"^ r/(0) = 2.3 x 10-^ r]'{0) = 0, 
L = 1. The solid line is the classical evolution, the dashed line is the evolution from the one- 
loop corrected equation ([4.51). Figure 2(b): One-loop contribution including the coupling g 
for the values of the parameters used in Figure 2(a). 

Figure 2(c): Velocity same values and conventions as in figure 2(a). 

Figure 2(d): {qmax{T)Y vs r for the same values as in figure 2(a). 

Figure 3(a): 77 vs r with g = 10"'' r]{0) = 2.27 x 10~^, r]'{0) = 0, L = 1. Same conventions 
for solid and dashed lines as in figure 2(a). 

Figure 3(b): One-loop contribution to the equations of motion including the coupling g 
for the values of the parameters used in Figure 3(a). 

Figure 3(c): Velocity |p, same values and conventions as in figure 3(a). 

Figure 3(d): (gmax(T))^ vs r for the same values as in figure 3(a). 

Figure 4(a): vs r with g = 10^^ r/(0) = 2.258 x 10~^ r]'{0) = 0, L = 1. Same 
conventions for solid and dashed lines as in figure 2(a). 

Figure 4(b): One-loop contribution including the coupling g for the values of the param- 
eters used in Figure 4(a). 

Figure 4(c): Velocity same values and conventions as in figure 4(a). 

Figure 4(d): {qmaxi'T))^ vs r for the same values as in figure 4(a). 
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